Fragmentation in Jets: Cone and Threshold Effects 
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We study the fragmentation of (light) quarks and gluons to hadrons inside a jet of cone size 
R. This allows for a more exclusive analysis of fragmentation than is currently the case. The 
shape of semi-inclusive cross sections in the hadron energy fraction z is described by fragmenting 
jet functions (FJFs), which we calculate in terms of R and the jet energy E. We introduce a new 
joint resummation to sum the double logarithms of R and 1 — z m the FJFs, which has a similar 
application to initial-state radiation at hadron colliders. Our results at next-to-leading logarithmic 
order indicate that the resummation of the threshold logarithms of 1 — z is already important for 
z > 0.5 and improves the convergence of perturbation theory. Our framework may be used to 
analyze LHC and RHIC data and to test and tune Monte Carlo event generators. 



Introduction. High-energy processes with an observed 
hadron in the final state can be described by factori- 
zing the short-distance (partonic) physics, which is per- 
turbatively calculable in QCD, from (universal) non- 
perturbative contributions, see e.g. Ref. In this con- 
text, the (unpolarized) fragmentation functions (FFs) 
D^{z,ii) [2] encode the information on the transition 
from an energetic parton i — {g, u,u, d, . . .} to a hadron 
h, which carries the fraction z of its energy, plus a re- 
mainder X. The knowledge of both perturbative and 
non-perturbative ingredients in factorization theorems is 
crucial to have control on theoretical predictions. For 
example, a better determination of the 6-quark FF was 
important for resolving a discrepancy between CDF data 
and theory for the pr spectrum of J/tp [3]. Furthermore, 
additional understanding of hadron production at high 
Pt in pp collisions is also required to determine more ac- 
curately the relative suppression of hadron spectra (jet 
quenching) seen in heavy-ion collisions |4j . 

Parameterizations for the FFs have been constrained 
by fitting to data for single-inclusive charged hadron 
production in e~^e~ at next-to-leading order (NLO) in 
perturbation theory [SI [3] . More recently, global analy- 
ses have been performed that also include semi-inclusive 
deep-inelastic scattering and/or pp, pp data from HERA, 
RHIC and the Tevatron [7J . To illustrate the current level 
of precision, the dominant D^^{z, /i — mz) is determined 
with uncertainties at the 10% level for z > 0.5 [6j. The 
FFs of the gluon and the non-valence quarks are known 
even less accurately. 

In contrast to the inclusive analyses listed above, the 
Belle collaboration is studying light-quark fragmentation 
in their on-resonance data, using a cut on the thrust event 
shape to remove the large ?>-quark background [S]. Here, 
the fragmentation takes place inside a (hemisphere) jet of 
invariant mass s, described by the fragmenting jet func- 
tions (FJFs) Qi{s,z,ii) that we introduced in Ref. [9]. 
Correspondingly, z should not be too small, to avoid con- 



tributions from hadrons outside the jets. For s 3> ^qcu^ 
the FJFs can be perturbatively matched onto the FFs [S] , 
i.e. — Jij (g) D'- where the convolution is in the 
momentum fraction z. In performing this OPE, we as- 
sume that parton and hadron masses ^ are negligible 
compared to the jet mass. The Wilson coefficients Jij 
describe the emissions from the parent parton at larger 
virtualities building up the jet while encodes physics 
at lower scales where hadronization effects are important. 

The one-loop Jij{s, z, were determined in Refs. [TUl 
[11], and allowed us to calculate in Ref. [IT] the cross 
section for e+e~ hX with a cut on thrust up to next- 
to-next-to-leading logarithmic order in Soft-CoUinear Ef- 
fective Theory (SCET)[T2]. There we found correlations 
between the thrust cut and z that are crucial for the ana- 
lysis of the Belle data. Of course this analysis is subject 
to non-perturbative power corrections, which can be siz- 
able at the Belle energy, and a comparison with data is 
desirable to assess their importance. 

Here we study fragmentation inside cone jets, defined 
with a cone jet algorithm. For simplicity we discuss the 
case of e"'"e~ — N jets, but our method can also be 
applied to pp collisions ^. Following Refs. [TSUlSj . the 
cone size is denoted by R and a cutoff A is applied on 
the energy in the region between jets ^. The jets are 
required to be energetic and well-separated: 



ta.T?{R/2), 



A 



tan2(i?/2) 
tan2(?/'/2)' 



< 1 



(1) 



^ The corrections due to the hadron mass are of order 

m^/(2:^Q^), where Q is the scale of the hard interaction. For 
light hadrons they are negligible unless z is small. 

^ For pp collisions the jet algorithms are defined in terms of <f> 
and the pseudo-rapidity ri and thus not azimuthally symmetric 
around the jet axis, introducing an additional hurdle. We assume 
the cancellation of Glauber gluons in factorization theorems. 

^ A may also be defined as a jet energy veto, which is equivalent 
up to 0{as). 
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where ijj is the minimum angular separation between jets 
and -Emin the minimum jet energy. Additional power 
corrections due to the jet algorithm are suppressed if R 
is not too small [TB], see also Eq. ([S]). SCET can then 
be applied to separate hard, collinear and soft dynam- 
ics, leading to factorization formulae as in Refs. [T51 114j . 
Schematically, at leading power. 
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(E, R)^j d$^tr[i/Ar5Ar] (S, R, z, ^x)\{Jt. 



(2) 



The hard function Hm describes the hard collision and 
the soft function Sm the soft radiation (both are matrices 
in color space). For each of the jets there is a jet function 
Ji describing the final-state collinear radiation, which is 
replaced by a FJF for the jet in which the hadron is ob- 
served [S]. The phase-space is denoted by d$Ar. The 
dependence on the renormalization scale ji cancels in the 
cross section up to the order that one is working and is 
used to estimate the uncertainty from higher-order cor- 
rections. Large logarithms of ratios of mass scales are 
summed by evaluating each of these functions at their 
natural scale, where they contain no large logarithms, 
and by running them to a common scale using their 
renormalization group equations. The FJFs with a cone 
restriction depend on E and R rather than s, and in 
the next section we calculate their matching onto FFs at 
NLO, 



A2 

4E?2tan2(i?/2) 



\i + oi 



(3) 



To avoid large non-perturbative corrections, R should 
thus not be too small ^. Since the other ingredients of 
the factorization theorem in Eq. ^ do not affect the 
fragmentation variable z, the shape of the cross section 
in z is completely determined by the cone FJF. (This 
is not the case when invariant masses are measured, as 
in Ref . [TT] , because invariant masses also receive a con- 
tribution from soft radiation.) The soft function is sen- 
sitive to two different scales [H] leading to non-global 
logarithms [17] but this only affects the normalization of 
the cross section. Furthermore, we would like to mention 
that the generalization of our framework to transverse- 
momentum-dependent and polarized FFs, can be applied 
as a tool to obtain information on the proton spin struc- 
ture by studying azimuthal asymmetries for hadron dis- 
tributions inside a high-p^ jet in transversely polarized 
pp collisions [T5] . 



* From Ref. jlBl non-perturbative effects are expected to scale 
like 1/R for small -R, which is consistent with the leading non- 
perturbative correction to the soft function. 



Calculation. We now calculate the one-loop matching 
coefBcients for the cone FJFs g^{R,z,ii) onto FFs. At 
one-loop the cone restriction is equivalent to 



s < mini : 



1 - z 



VI 



AE^ tan2(i?/2). 



(4) 



where s is the invariant mass of the jet, z the fragmenta- 
tion variable, E the jet energy and R the cone size. [Note 
that E and R appear in the combination E^ta.n^{R/2) 
which is invariant under boosts along the jet axis.] We 
may therefore obtain the matching coefficients from the 
bare results for the standard FJFs in Ref. Hj. We find, 
using the MS scheme, 



Jqg{E,R,Z,fl) 

2(27r)3 
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+ Pgg{z)L + Jgg{z) 



S{l-z) + 



Pggiz)L + Jggiz)] 

UsCa 



*(l-=)(i"-y 



+ Pgg{z)L + Jgg{z) 



Jgq{E,R,z,^) _ a^Tp 
2(27r)3 ~ TT 



Pgg{z)L + Jgg{z)\ , 



(5) 



with the splitting functions in the convention of Eq. (3.7) 
in Ref. T^. Anti-quarks have the same coefficients as 
quarks, and Jqq and Jqqi only start at two-loop order. 
In Eq. ^ 



L = In 



2£;tan(i?/2) 



(6) 



and 



1 



JM = -{l-z) + 
'^ggi'^) ^ 2(i-z+ 

Jgqiz) ^ Zil - Z) + Pqgiz) 



Pqq{z) In 2 
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ln(l - z 



z < 
) z> 
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ln(l - z) z > i 



(7) 



We have cross-checked these results with expressions 
for the jet functions Ji with a cone restriction in Ref. [T4] 
(there called "unmeasured jet functions"). The FJFs Gi 
and jet functions Ji have the same anomalous dimension 
[5], consistent with Eq. ([2]), which provides a check on 
the UV-divergent terms that we encounter in deriving 
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FIG. 1: The FJF for u -)■ 7r+ at LL and NLL order for E = 100 GeV and R = 0.4. The bands show the perturbative 
uncertainties from varying the scale in the matching onto DJ up and down by factors of two, and then evolving QZ to 
fj, = 2E tan(7?/2). For reference we include the LO result (where fj, is the jet energy E) in the left panel. In the middle panel 
the same curves and bands relative to this LO result are shown. The dotted lines correspond to the central values without 
threshold resummation, for which the uncertainties are shown in the right panel. 



Eq. ([5]). The finite terms were checked with the momen- 
tum sum rule HI [TT] 



"'0 



(8) 



2{2nf Jq[R,y). 



and a new sum rule 
»i 

dz [Jqq{R, Z, - JqqiR, Z, ^)] 

(9) 

This new relation follows from quark-number conserva- 
tion in the perturbative calculation of the FJF tj*. 

Threshold Logarithms. Closer inspection reveals that 
Jqq and Jgg contain large threshold logarithms. For ex- 
ample, 



Jqq{E, R, z, fl) _ UsCf 



2{2nf 



2L 



2{ 



ln(l - z) 



24. 
0(1 



, (10) 



V I- z 

which by Eq. ^ leads to the large (double) logarithms 



g^{E,R,z,fi) 
2(27r)3 



= ln^ 



2(1 - z)£;tanf 



D'^{z) + ... (11) 



We may thus sum both the logarithms of R and the 
threshold logarithms by evaluating the FJF at /i = 
2(1 - z)£'tan(i?/2). (In the plots we evolve the FJF to 
/i = 2i? tan(_R/2), to remove the z dependence from this 
scale.) Fig. [1] shows the improved convergence of pertur- 
bation theory arising from threshold resummation, which 
will be discussed more extensively in the next section. We 
stress that we do not need to perform a threshold expan- 
sion; therefore our results are also valid away from the 
threshold region. 

This improvement holds to all orders in perturbation 
theory, since we find that the logarithms of R and 1— z are 
tied together through the anomalous dimension of jT^g ^ . 



Using Eq. ^ this can be obtained from the anomalous 
dimension of the (unmeasured) jet function [T31 [Tlj and 
the FF m. 



TiT qq 



-2r„ 



LS{1 - z) 



(1-^)h 
+ jjS{l-z)+0{l-z). 



(12) 



Here F, = asCp/TT + 0{al) is the (quark) cusp anoma- 
lous dimension. The non-cusp anomalous dimension jj 
only starts at two loops, which is refiected in the ab- 
sence of single logarithms in Eqs. (10) and (111. In mo- 



ment space one can see explicitly how the logarithm of 
R and the 1/(1 — z)^ combine, since the iVth moment of 
1/(1 -z)+ is - \nN--/E- 

The above discussion applies to J^gg as well. By con- 
trast, the off-diagonal J^qg and J7gg do not contain thresh- 
old logarithms and their contribution should thus be eval- 
uated at the appropriate scale /i = 2iJ tan(i?/2). 

Our new joint resummation can also be used in the 
standard FJFs Qf{s,z,^) of Refs. [5HTT]. Here the scale 
choice fx^ = (1 — z)s sums both the logarithms of s and 
1 — z. We have checked that this significantly improves 
the convergence of perturbation theory at large z. 

Similarly, for initial-state jets described by beam func- 
tions Bi{t, X, fi) IH], we find that the logarithms of t and 
1 — X can be simultaneously resummed by evaluating 
the flavor-diagonal terms for matching beam functions 
onto PDFs at the scale //^ = (1 - x)t. At the LHC, the 
momentum fraction x is typically small, suggesting that 
threshold resummation is less important than in the case 
studied here. However, it has been argued that the log- 
arithms of 1 — a: are dynamically enhanced through the 
shape of the PDFs [22] . A detailed study is left for future 
work. 



^ This requires the all-orders structure of the jet function anoma- 



lous dimension, which is known for hemisphere invariant mass 
jets, but not established for jet algorithms. 
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FIG. 2: The FJF for u,d,g ^ n+ is shown at NLL order for E = 100 GeV. The curves are for R = 0.2,0.4,0.6,0.8 and are 
shown relative to the R = 1 result. 



Results. We now show some numerical results for 
fragmentation inside a cone jet. Our order counting is 
shown in Table |lj As input we use the HKNS FFs _6; at 
LO or NLO order and their value for a^. 

We first consider a fixed cone R = 0.4 and study 
the convergence of our resummed result. Fig. [T] shows 
the central values and perturbative uncertainties for the 
M — > Tr"*" FJF at LL and NLL order. In the middle 
and right panels these results relative to the LO curve 
are shown, making the separation between the central 
curves and the size of the uncertainties more visible. As 
expected, the cone restriction suppresses the small z re- 
gion because at large z the hadron carries off most of the 
energy and there is not much additional radiation that 
would be subject to the cone restriction. As is clear from 
comparing the dotted with the solid curves, the effect of 
threshold resummation becomes important ad, z ^ 0.5. 
Including threshold resummation improves the conver- 
gence and leads to smaller uncertainties. At small z the 
difference between LL and NLL is not fully captured by 
the uncertainty band. However, in this region the uncer- 
tainties on the FFs (which are not shown here) dominate. 

We have studied similar plots for d — tt"*" and g — >■ tt'^ . 
Here the convergence is not as good when the dominant 
mixing contribution involving (and D'^ ) turns on. 

For (7 — 7r+ and d — >■ 7r+ this happens at NLO and 
NNLO, respectively. 

In Fig.[2]we study the effect of varying the cone size R. 
We stress that the normalization of the curves for differ- 
ent R should not be compared because other ingredients 
in the factorization theorem affect the normalization of 
the cross section (the shape in z is entirely determined by 
Q). For u — > 7r+ and d — > 7r+ the small z region gets sup- 
pressed relative to the large z region as one reduces the 





Jij 7 


r 13 


LO 


0-loop 


1-loop 


NLO 


1-loop 


2-loop 


LL 


0-loop 


1-loop 1-loop 


NLL 


1-loop 1-loop 2-loop 2-loop 



cone size, as expected. For g — tt''", the bump-shaped 
enhancement at intermediate z is mainly due to the J^gq 
contribution. 

Conclusions. Our framework leads to a reliable de- 
scription of fragmentation within an identified jet, pro- 
viding the tools for a more exclusive study of fragmen- 
tation. This can be used to reduce the background 
from other processes, to help identify underlying partonic 
structures and enables novel tests of the universality of 
FFs. Our setup provides accurate analytical predictions, 
where large logarithms are properly summed, for single 
hadron spectra in a jet which can be used to tune Monte 
Carlo event generators. This will become more and more 
relevant as the determination of the fragmentation func- 
tions becomes more accurate. We have for the first time 
established the connection between the resummation of 
logarithms for exclusive jet production and threshold re- 
summation, and have shown the numerical importance 
of their interplay in the case of jet fragmentation. 
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TABLE I: Order counting for the matching coefficients, the 
non-cusp and cusp anomalous dimension, and the as running. 
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